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Abstract 

We extend previous work on the arrival time problem in quantum mechanics, in the framework 
of decoherent histories, to the case of a particle coupled to an environment. The usual arrival 
time probabilities are related to the probability current, so we explore the properties of the current 
for general open systems that can be written in terms of a master equation of Lindblad form. 
We specialise to the case of quantum Brownian motion, and show that after a time of order the 
localisation time the current becomes positive. We show that the arrival time probabilities can 
then be written in terms of a POVM, which we compute. We perform a decoherent histories 
analysis including the effects of the environment and show that time of arrival probabilities are 
decoherent for a generic state after a time much greater than the localisation time, but that 
there is a fundamental limitation on the accuracy, 6t, with which they can be specified which 
obeys E6t >> h. We confirm that the arrival time probabilities computed in this way agree with 
those computed via the current, provided there is decoherence. We thus find that the decoherent 
histories formulation of quantum mechanics provides a consistent explanation for the emergence of 
the probability current as the classical arrival time distribution, and a systematic rule for deciding 
when probabilities may be assigned. 
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I. INTRODUCTION 



Questions involving time in quantum theory have a rich history, and there is still much 
debate about their status [1]. Quantities such as arrival and dwell times, despite being 
measureable [2] , still lack concrete grounding within "standard" quantum theory, and there 
has been considerable interest in understanding these quantities in the framework of various 
interpretations of quantum theory. Arrival times, in particular, have attracted much interest, 
as the natural procedure of quantising the appropriate classical quantity gives rise to an 
operator which is not self-adjoint and thus, in standard quantum theory at least, cannot 
easily be considered as an observable. 

Despite these difficulties, if one considers a free particle in an initial state p = {ip) (V'l 
localized in a; > consisting entirely of negative momenta, and asks for the probability 
p(ti,t2) that the particle crosses the origin during the time interval [^1,^2] there is a semi- 
classical answer, given by the difference between the probability of being in x > at the 
initial and final time [I]. Defining P{t) = 6{x{t)), 

p{h,t2) = Tr(P(ti)p) - Tr(P(t2)p) = - / ' dtTT{Ppt) 

dtJ{t), (1.1) 



where 



is the standard quantum mechanical probability current. (We denote the state at time t by 
Pt-) This can also be written in the following two forms, 

p{tut2) = r dt [ dpdq^-^^^Wt{p,q), (1.2) 

where Wt{p,q) is the Wigner function, defined later, and 

p{t,,t2) = Tr(Cp) 

C = r dt^^ {pS{x{t)) + 6{x{t))p) 



tl 



2m 

P{h)-P{h). (1.3) 



These expressions agree with the classical result, with the Wigner function W replaced by 
the classical distribution function w, provided that the classical trajectories are straight lines. 



They are also correctly normalized to 1 as ^2 — C)0 with ti = 0, but they are not generally 
positive, even for states consisting entirely of negative momenta. This genuinely quantum 
phenomenon is called backfiow and arises because the operator C, positive classically, has 



negative eigenvalues pHH]. This means we cannot generally regard Eq.(l.l) as providing 
an acceptable arrival time distribution. There is, in addition, a more fundamental problem 



with Eq.(l.l), which is that probabilities in quantum theory should be expressible in the 
form |6] 

via) = Tr(P,p). (1.4) 
Here is a projection operator, or more generally a POVM, associated with the outcome 



a. Eq.(l.l) cannot be expressed in this form, and we therefore conclude that it is not a 
fundamental expression in quantum theory, so must be the result of some approximation. 



An interesting clue as to how to improve on Eq. (1.1) is provided by the expression for the 



current in terms of the Wigner function, Eq.(1.2). Negativity of the Wigner function is a 
neccessary condition for the negativity of the current (in the sense of backfiow, as discussed 
in [7]), but it is known that evolution in the presence of an environment typically renders the 



Wigner function positive after a short time [8]. This suggests that something like Eq.(l.l) 
may be an acceptable, if still heuristic, arrival time distribution for a system coupled to an 
environment, at least after some time. The first aim of this paper is to derive the correct 



analogue of Eq.( 1.1 ) for a system coupled to an environment, and to prove that it can indeed 
be regarded as an arrival time distribution after a short time. In particular we show that the 
arrival time probabilities computed in this way may be written as the trace of the density 
matrix times a POVM. 

However although this may allow us to interpret the current as the arrival time distri- 



bution, we are still left with the task of deriving Eq.(l.l) from some more fundamental 
quantity. This has previously been achieved, in some approximation, in Refs.[71 19] in the 
context of the decoherent histories approach to quantum theory [IDHIll, for the case of a 
free particle, and the derivation was shown to hold for states and intervals exhibiting suffi- 
cient decoherence. In decoherent histories the probability of arriving in an interval [^1,^2] is 
approximately given by 

pih,h)=TT{CpC^), (1.5) 



with C given by Eq.(1.3), and this reduces to Eq.(l.l) under the somewhat special condition 
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of decoherence. However it was also shown that there exist states for which decoherence of 
histories cannot be obtained for reasonable levels of coarse graining, and for which therefore 
we cannot assign a time of arrival probability distribution. States exhibiting backflow are 
good examples of this. 

Lack of decoherence for a general state is a situation frequently encountered in the liter- 
ature on decoherent histories. The solution lies in the observation that realistic systems are 
always coupled to an environment, and that as such they are most fundamentally described 
by open system dynamics. Because by definition an environment consists of degrees of free- 
dom about which we have no knowledge, and over which we have no control, it is natural 
to coarse-grain over these degrees of freedom. Such coarse-graining generally results in the 
recovery of approximate classical behaviour, and thus decoherence. We therefore anticipate 
that decoherence of histories corresponding to arrival times may be achieved for a generic 
state, provided that particle is coupled to a suitable environment. The second aim of this 
paper is to examine this scheme. We show that the probabilities calculated in this way are 
approximately decoherent, and approximately equal to those computed via the analogue of 



Eq.(l.l), valid in the case of an environment. 

Once we introduce an environment, however, the correspondence is with a classical 
stochastic theory, rather than with a deterministic one. The classical trajectories may now 
cross the origin many times due to fluctuations. The arrival time distribution for such a 
classical theory was computed in Ref . [13] , and is given by 

Pih,t2) = r dt j dp j dq^—^6{q)wl{p,q), (1.6) 



where q) is the initial phase space distribution, evolved with a type of restricted prop- 
agator valid in the case of an environment, and with boundary conditions, [15] 



^/;[(p,0) = 0, for p>0. (1.7) 

In Ref.|15] Halliwell and Zaflris presented a decoherent histories analysis of the cor- 
responding quantum case. Although the conclusion they reach is sensible their analysis 
actually contains a small error. This was the result of a lack of appreciation of the role of 
the quantum Zeno effect in these calculations. We will show in this paper how the analysis 
can be corrected to give a correct expression for the arrival time probability in decoherent 
histories. We shall be interested in a particular limit of this general case where the stochastic 



trajectories are sharply peaked about the deterministic trajectories we would have in the 
absence of any environment. In this case the restricted propagator may be replaced with an 
unrestricted one, and we will show that quantisation in this limit yields an expression of the 



form Eq.(1.2), but with the Wigner function evolved in the presence of an environment. 

It will eventually be necessary to specialise to a particular model of system-environment 
coupling in order to obtain quantitative results, but we will begin by considering general 
models which can be written in terms of a master equation of Lindblad form. This kind 
of system has been extensively studied in the decoherence literature, and the properties 
of such evolution, such as suppression of interference effects, loss of entanglement and an 
approximate recovery of classical behaviour have been discussed in [8], ITGHTS] . 

This paper is arranged as follows, in Section [TT] we explore some properties of the arrival 
time distribution for general open systems of Lindblad form, our aim being to derive the 
corrections to the current resulting from the environmental dynamics. In Section III we 



discuss quantum Brownian motion and specialise the results of Section [IT] to this case. In 



Section IV we derive some properties of the arrival time distribution for quantum Brownian 
motion, and in particular we prove that current becomes positive after a finite time. In 



Sections \V\ and [VI] we discuss the decoherent histories approach to quantum theory, and the 

we 



expressions for the arrival time distribution that arise in this context. In Section VII 



then examine whether histories corresponding to arriving in different intervals of time are 
decoherent. We summarise our results in Section IVIIII 



II. ARRIVAL TIME FOR OPEN QUANTUM SYSTEMS 

We consider an open quantum system consisting of a free particle coupled to an environ- 
ment with master equation of the following Lindblad form |19j . 

I = P] - ^ E ({^1^- P} - 2^.P^]) ' (2-1) 

j 

where Lj are the Lindblad operators, and H = Hq+Hiis the free hamiltonian plus a possible 
extra term that depends on the Lj. (This term gives rise to frequency renormalisation in 
oscillating systems.) Specific forms for the Lj may be computed for particular models, but 
for the moment we leave this general. The Lindblad equation represents the most general 
master equation for a Markovian system which preserves the properties of the density matrix, 
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in particular positivity. 



We now extend the derivation of Eq.(l.l) to this system. The probabihty of crossing 
during the interval [ti,t2] is 



p{h,t2) = Tr(Ppi J - Tr(Pp 



t2y 



dm{Ppt 



ti 



-1 

2m 



t2 



dtTr {\pS{x) + S{x)p]pt) — 



t2 



dtTr{[H,,P]pt 



+ 1 E ^^Tr {[L], P]L,p, + L][P, L,]p, 



(2.2) 



where P = 6{x). The first term is the standard current expression, although with the state 



evolved according to Eq.(2.1), and we therefore recover Eq.(l.l) when all the Lj are 0. The 
second and third terms depend on the Lindblad operators, Lj and thus on the form of the 
system-environment coupling. 

To proceed further we specialise to the case where L is a linear combination of x and p, 
L = ax + ibp, where a and b are real constants. The master equation Eq.(2.1) is then 

.2 



dpi 
dt 



1 ^2 
--[Ho,pt] -iab[x, {pt,p}] - yl^, [x,pt]] - — [p, [p,pt]]. 



(2.3) 



Note that this equation is also identical in form to the exact master equation for a particle 
in a gas environment, given in j20] . 

To derive the arrival time distribution we could simply substitute L = ax + ibp into 



Eq.(2.2), but the algebra is somewhat clumsy, and there are a number of terms which 
vanish for reasons not immediately apparent from this expression. An equivalent approach 



is to start from Eq.(2.3), multiplying this expression by P and taking the trace. Because the 
second and third terms on the right of this expression have the form [x, A] their contribution 
is proportional to 

Tr([x, i]P) = Tr(i[P,x]) = 0,, 



so that only the first and last terms on the right hand side of Eq.(2.3) contribute. The final 



term on the right of Eq.(2.3) gives a contribution 
52 u2 



-Tt{[p, [p,pt]P) 



Tr([p, [P,p]]pt) = -thb'TT{[p6{x) - 6{x)p]pt), (2.4) 
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so we arrive at the expression, 



p{UM) = ^ dm {[p6{x) + 6{x)p]pt) 



rt2 

-ihh^ I dm {\p5{x) - 5{x)p\pt) . (2.5) 

The second term in this expression has a somewhat unusual form and its significance is 
not immediately clear. We will see below that this term is related to diffusion in position. 
There is a connection here to a recent paper by Genkin, Ferro and Lindroth [21], in which 
the authors sought to examine the effects of an environment on the arrival time distribution. 
In that paper the authors used the standard expression for the current valid in the unitary 



case, Eq.(l.l), ignoring the extra terms that arise because of the environment. Although this 
may be a good approximation when we can neglect the effects of dissipation, it is clear that 
there may be significant corrections to the current for strong dissipation. They also pointed 
out that these corrections are equivalent to the presence of extra terms in the continuity 
equation, although they did not compute these explicitly. For the sake of completeness, and 



also because it helps to understand the nature of the extra terms in Eq.(2.5) we will derive 
them here. 



To derive the continuity equation we multiply Eq.(2.3) by 5{x — x) and perform the trace. 
If we neglect the final three terms on the right we arrive at the standard result, 

^(x,x) + |^(a;,t) = 0, (2.6) 



with J{x,t) defined by an obvious extension of Eq.(1.3). Turning to the extra terms that 
result from the inclusion of the environment, the second and third terms vanish in exactly 
the same way as for the current above, and the correction term is therefore given by, 

-^TY([p,[p,p,Mx-x))=2;iV^(x,x), (2.7) 

so the continuity equation now reads, 

^{x,x) + ^ (J(x,t) + JD{x,t)) = 0, (2.8) 

where we have identified the diffusive current. 



Joix^t) = -2h%'^^ix,x). (2.9) 
ox 
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This is a specific example of a modification to a conservation equation resulting from open 
system dynamics. For a more general discussion of these issues see [22] • 

The analysis presented above can also be phrased in terms of phase space distributions. 
The Wigner function corresponding to pt is defined as [23] 



(2.10) 



The Wigner transform of the master equation Eq.(|2.3|) is 

dWt 



p dWt 



+ 2h^ah 



+ 



+ 



(2.11) 



dt m dq dp 2 dp"^ 2 dq^ 

The first term is the standard unitary term, whilst the second term represents dissipation, 
and the third and fourth terms represent diffusion. The arrival time distribution can be 



written in terms of the Wigner function by taking the Wigner transform of Eq.(2.5) 



p{ti,t2 



ti 



dt / dpdq 



i-n) lYlf dWf 

^-^mWt{p,q) + V^(g)^(p,g) . 

m 2 oq 



(2.12) 



Eqs.(2.5) and (2.12) are the sought for generalisation of Eq.(1.2) for the case of a particle 
coupled to an environment. 



III. QUANTUM BROWNIAN MOTION 

Quantum Brownian motion [161 113 El] is a commonly used form for the environment of 
an open system, partly because it is exactly solvable, and partly because in many cases it is 
a good approximation to a realistic environment. For the quantum Brownian motion model 



[T7] we have Eq.(2.3) with 



1. 

Here D = 2m'ykT, where T is the temperature of the environment, and 7 is a phenomeno- 



j2D/h^, h=^=. (3.1) 



logical damping constant, see [21]. Eq.(2.3) may be written as 



Although we could in principle work with this general case, it is useful to specialise to 
the case of negligible dissipation. There are two reasons for this, the first is that the analysis 
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is considerably simplified, and this helps us to see the important effects more clearly. The 
second is that we have a particular aim in mind here, and that is to understand how a sensible 
classical result emerges from the quantum case. The classical case we have in mind is that 
of a heavy particle following an essentially classical, deterministic trajectory, but subject to 
small quantum fluctuations. We therefore restrict our analysis to timescales much shorter 
than the relaxation time 7"^, and so we can drop the final two terms in the master equation 
above. 

This master equation may then be solved in terms of the propagator 



Pt{x,y) = J dxodyoJ{x,y,t\xo,yo,0)po{xo,yo) (3.3) 

/ in \ ( tin 
J{x,y,t\xQ,yo,0) = j I ^ ~ ^0)^ ~ ~ ?/o)^] 



Taking the same limit in the equation for the Wigner function, Eq.(2.11), gives 



dt m dq dp^ 
Evolution of the Wigner function may also expresed in terms of a propagator |16] . 



Wt{p,q) = / dpodqoK{q,p,t\qo,po,0)Wo{po,qo) (3.6) 



2 o I ... Pot 



K{q,p,t\qo,Po,0) = N exp \ -a{p - po) - /3 { q - qo 



m 



2 



+e(p-Po) ( g-go- — ) ) , (3.7) 



m 



where A^, a, (3, and e are given by 

1 ^ 3m^ 3m / 3m^ ^ ^^'^ 



a = -—. 
Dt 



^ 3m^ 3m / 3m^ \ ^ 



For later convenience we note that we can make the simple change of variables here, 

Qo ^ Qo — Pot/^, so that 

Wtip,q) = j dpodqoKiq,p,t\qo,po,0)Wo{po,qo) (3.9) 
K{q,p,t\qo,po,0) = N exp {-a{p - po)"^ - l3 {q - qof + e{p - po){q - qo)) 

W{p,q) = W{p,q-pt/m). (3.10) 



This form makes it clear that the evolution consists of two effects. The first is a shifting 
along the classical trajectories, whilst the second is a spreading in phase space. 

It is useful to consider this process in more detail. In the presence of an environment the 
width of the momentum distribution becomes time dependent, and we have, 

(Ap)? = {Ap)l + Dt, 

where (Ap)o is the momentum width of the initial state. We recognise an immediate difficulty 
here. Even for initial distributions consisting entirely of left moving momenta, Wt will 
develop support on p > under evolution. This means we cannot strictly regard the current 
as an arrival time distribution, since a typical trajectory will now cross the origin many 
times. Differently put, the arrival time distribution is strictly a measurement of the first 
passage time, and this is no longer equal to the current. (Note that although this result is 
similar to the problems created by backfiow, the reasons behind it are very different. The 
spreading of momentum induced by evolution in an environment is a purely classical effect.) 

However, all is not lost. Although the current is no longer strictly the arrival time 
distribution, it may be a very good approximation to it. This is because the deviation 
of the current from the "true" arrival time distribution will be related to the probability 
that the state has the "wrong" sign momenta. This means that, provided we are in the 
"near-deterministic" limit, 

{Apf,«pl (3.11) 

where po is the momentum around which the initial state is tightly peaked, the current will 
still be a very good approximation to the true arrival time distribution. The timescale on 
which this analysis breaks down is given by the "stochastic" time, 

rs=pl/D. 

After this time we must therefore revert to using the exact expression for the arrival time 



given by Eq.(1.6). We see from the definition of D that. 



for the states we are interested in. This means the stochastic time is much longer than 
the relaxation time 7"^, so that working in the near-deterministic limit does not impose any 
further constraints compared with neglecting dissipation. 
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Returning to our arrival probability, in the limit of negligible dissipation Eq.(2.5), be- 
comes, 

p{t2,ti) = ['dtJ{t) (3.13) 

J{t) = f dqdp^—^5{q)Wt{p,q). (3.14) 
J m 



This expression is now identical to the unitary case, Eq.(l.l), but with the Wigner function 
evolved under quantum Brownian motion. 

We now turn to the question of whether inclusion of an environment helps ensure the 



positivity of Eq.(3.13). 



IV. PROPERTIES OF THE ARRIVAL TIME DISTRIBUTION IN QUANTUM 
BROWNIAN MOTION 

In the introduction we noted that evolution in an environment typically renders Wt 
positive after a short time. We now wish to examine the effect of this on our candidate 



arrival time distribution Eq.(3.13). To this end we introduce the notation [8j 

(4.1) 




and also the class of Gaussian phase space functions 

where A is a 2 x 2 positive definite matrix, with determinant \A\. g{z; A) will be a Wigner 
function if and only if 

1^1 > J- (4.3) 



A useful result is that 



j d^zg{zi - z; A)g{z - Zs; B) = g{zi - Zg; A + B). (4.4) 



In this notation we can write the propagation of the Wigner function, Eq.(3.6), as 



Wt{z) = j dz'g{z - z'- A)W,{z'), (4.5) 
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where 



A = Dt 



2 t/m 
t/m 2tV3m2 



Since 



after a time 



1^1 



DH 



2+4 



16 



1/4 



3m2 
2mh^ ^^"^ 



D 



3 

16 



1/4 



n 



(4.6) 



(4.7) 



(4.8) 



g{z — z'; y4) will be a Wigner function, and thus Wt{z) will be positive because it is equal to 
the convolution of two Wigner functions. Here ti = ^J2mh/D is the localisation time. 



This is a useful result. Expressing the current Eq.(3.14) in this new notation, 



Jit) 



dz'-^6izi)Wtiz) 



m 



(4.9) 



we see that after the time given in Eq.(4.8), because Wt > the current Eq.(4.9) will be 
positive if the state consists purely of negative momenta. This means that after this time 



Eq.(3.13) will be a positive arrival time distribution. This holds provided times involved are 



much smaller than r^, as per the discussion below Eq.(3.10). 

Now we turn to examining the properties of the current. We wish to find an expression 
for the current in the form Tr (Pp), where P is a projector or POVM. This would allow 



us to express the heuristic arrival time distribution, Eq.(1.2), in the same form as other 



probabilities in quantum theory. Starting from Eq.(4.9) it is useful to write, 

-zo) 



Jit) 



dz / dz'^ ^5(zi)^(z-z';A)Wo(z') 



m 



(4.10) 



using Eq.(4.5). We can deconvolve the propagator into two gaussians using Eq.(4.4), in 
particular we will let A = Aq + where Aq is a minimum uncertainty gaussian, and B is 
the remainder. 



Aa = h 



■ 

1/4:/' 



B = Dt 



2 t/m 
t/m 2tV3m2 



-Ao. 



(4.11) 



(4.12) 
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Here s is some real number. Using the convolution property Eq.(4.4) we can write the 



current Eq.(4.9) as 



At) 



az 
dz" 



dz^^6iz^)giz-z";B) 



m 



dz'Wo{z')g{z" -z';Ao) 



dz^^6{z^ + Zot/m)g{z - z"; B) 



m 



(4.13) 



where Q{z) is the Q-function |25], and we have undone the change of variables implied in 



Eq.(3.10). The Q-function can be written as 

g(z) = - (z|p|z) . 



TT 



We can therefore express the current as, 



J{t) = Tr(Fp), 



(4.14) 



where 



^ j dz" \z") {z"\ (^j dz^^-^6{zi + Zot/m)g{z - z"; B) 

( — Zq) 

dzP^ 6(zi + ZQt/m), 

m 



(4.15) 



and we have defined the POVM 



Pz = ^ j dz" \z") {z"\ g{z - z"; 5), 



(4.16) 



which is clearly a phase space operator localised around z. F is therefore a smeared version 
of the object used to compute the current classically, ~p6{xt)/m. This holds for times 

o X 1/4 

- ] n<t « r,. (4.17) 
Assume for a moment that B = 0, and so Pz = |z) (z|. The current would then be given 



by 



At) 



dz- 



m 



6{zi - ZQt/m)Q{z). 



(4.18) 



Since the Q-function is positive by construction the current computed in this way will be 
positive to the extent that the Q-function has vanishing support on p > {Q{z) cannot 
strictly vanish for p > but it can be exponentially small, and this suffices here. See the 
comments below Eq.( 3.10) ).) For B > the Q-function is simply smeared further, and this 
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will preserve the property of positivity. In fact one could imagine postulating Eq.(4.18) as 
the definition of the arrival time even in the unitary case, since it clearly satisfies all the 
conditions required. 



The probability of arriving in the interval [^1,^2] is given by the integral of Eq.(4.14), and 
because Pz is time dependent (via B) this will not in general have a simple form. However 



there is a separation of time scales here, the time scale of evolution of Pz is seen from Eq.(4.7) 
to be Ti. So, if t2 — ^1 << Ti we have approximately 



where 



dtTT{Fp) ^ / rfzTr 



t2 



dt- 



m 



6{zi + ZQt/m) 



Tr(Ep), (4.19) 



E = j dzPz[eizi + zoti/m) -9{zi + zot2/m)] (4.20) 

is a POVM representing arrival at x = between ti and t2- Note the close similarity between 
this object and the operator Eq.(1.3), but that crucially this phase space operator is positive 
for zo < (ie p < 0). In this case therefore, the effect of the environment is simply to smear 
the unitary result over a region of phase space, given by computed at ti. 



Eqs.(4.19) and (4.20) form the first significant result of this paper. Eq.(4.19) expresses 



the heuristic arrival time distribution, Eq.(1.2), as the trace of an operator times a POVM, 
and thus has the same form as standard expressions for probability in quantum theory. The 



POVM, Eq.(4.20), arises because the environment effectively measures the system. 



We have therefore discovered a range of times for which the current, Eq.(4.9), gives a 
positive arrival time distribution. After a time of order r/ interference effects vanish and 



we can regard Eq.(4.9) as the arrival time distribution, which can also be written in the 



form Eq.(4.14). Eventually, however, on a time scale of order the environment causes 
diffusion in momentum to such an extent that the trajectories of the particle are no longer 
sharply peaked around the classical trajectory computed in the absense of an environment. 
Since the particle is still behaving classically after this time there will exist an arrival time 



distribution of the form Eq.(1.6), but our simpler expression Eq.(4.9) will no longer be a 
good approximation to this. It is easy to show that Ts/ti = EtiH where E is the energy of 
the initial state, and thus these expressions are valid for an large interval if Eti » h. 
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V. THE DECOHERENT HISTORIES APPROACH TO THE ARRIVAL TIME 
PROBLEM 



A. The Decoherent Histories approach to Quantum Theory 

We begin by briefly reviewing the decoherent histories approach to quantum theory. More 
extensive discussions can be found in the references |T0HT1]. 

Alternatives at fixed moments of time in quantum theory are represented by a set of 
projection operators {Pa}, satisfying the conditions 

E^'^ = 1 (5.1) 

a 

PaPb = SabPa, (5.2) 

where we take a to run over some finite range. In the decoherent histories approach to 
quantum theory histories are represented by class operators Ca which are time-ordered 
strings of projections 

Ca = PaM ■ ■ ■ PaAh) (5.3) 

or sums of such strings [26j. Here the projections are in the Heisenberg picture and a denotes 
the string (ai, ■ ■ ■ a„). All class operators satisfy the condition 

J2Ca = l. (5.4) 

a 

Probabilities are assigned to histories via the formula 

p{a) = Tr {CapCi) . (5.5) 

However probabilities assigned in this way do not necessarily obey the probability sum rules, 
because of quantum interference. We therefore introduce the decoherence functional 

D(a,/3) = Tr(c„pCj), (5.6) 

which may be thought of as a measure of interference between pairs of histories. We require 
that sets of histories satisfy the condition of decoherence, which is 

D{a,(3)=0, a ^(3. (5.7) 

This ensures that all probability sum rules are satisfied. 

15 



We note briefly that when there is decoherence Eq.(5.4) and Eq.(5.7) together imply that 
the probabihties p{a) are given by the simpler expressions 



=Tr(C„p). (5.8) 

Decoherence ensures that q{a) is real and positive, even though it is not in general. In this 
way decoherent histories may reproduce probabilities of the form Eq.(1.3). 

B. The Decoherent Histories aproach to the arrival time problem 

We turn now to the deflnition of the arrival time problem in decoherent histories. In 
Refs.[7l |9] the class operators corresponding to a flrst crossing of the origin in the time 
interval [tk-i,tk] were computed to be 

Ck = P{tk)P{tk-i) ■ ■ ■ P{t2)P{U), (5.9) 

for k > 2, with Ci = -P(ti) and where P{t) = 6{—Xt). These clearly describe histories which 
are in x > at times ti,t2, - ■ ■ tfc-i and in x < at time tk, so, approximately, describe a flrst 
crossing between tk-i and tk- We are prevented from taking the time between projections 
to zero by the Zeno effect . This point was extensively discussed in [3 E] • 

In Ref [9] these class operators were then simplifled with the help of the following semi- 
classical approximation 

P{t^) ■ ■ ■ P{t2)P{hm ^ P{tnM). (5.10) 

This is the statement that given that the state is in x > at t„, it must also have been x > 
at all previous times (we restrict attention here to states with p < 0). It is clear from the 
path integral representation of the propagator that this is only true semiclassically. From 
this we obtain the class operators for flrst crossing between tk-i and tk as 

Ck ~ Pitk)Pitk-i) (5.11) 
= P{tk-i) - P{tk)P{tk-i) 

^ P{tk-i) - P{tk), (5.12) 

where we use the semiclassical approximation again to arrive at the last line. In the case 



of a free particle without an environment, this class operator reproduces Eq.(l.l) under 
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the conditions of decoherence. This is because, assuming decoherence, the arrival time 
probabihty computed from decoherent histories is 

p{t2,h) = Tr(CpCt) = Tr(Cp), (5.13) 



using Eq.(5.8 ). 



VI. DECOHERENCE AND THE SEMI-CLASSICAL APPROXIMATION FOR 
ARRIVAL TIME 

In the previous section we introduced the arrival time problem in decoherent histories. 



The key step in making contact with the simple heuristic result of Eq.(l.l) was the ap 



proximation of the class operators Eq.(5.9) by Eq.(5.12). We now want to explore this 
approximation, and prove that it is indeed valid in some interesting limit. 

For a free particle, the error in this approximation comes from the finite width of the 
wave packet. This means that 

Pit,)\ij)^\^P), (6.1) 

even for times, tk, significantly before the classical crossing time. However it is clear that 
the largest error will come from the last projection that preceeds the classical crossing time. 
Since the spacing between projections has to be greater than the Zeno time, which is also 
the time scale on which a wavepacket may be said to "cross" the origin, only one or two 
projections stand any chance of significantly disturbing the wavepacket. 
Writing this in terms of the density matrix, 

TTiP{tk+i)P{h)pP{tk)) = Tr(P(t,,+i)p) - Tr(P(t,+i)P(4)p) - Tr(pP(tfc)P(tfc+i)) 

+Tr(P(4+i)P(tfc)pP(tfe)) 
^ Tr(P(tfc+i)p), 

where the last line defines the semi-classical approximation. Noting that, with tm > tk 

|Tr(P(OP(4)p)|' < TT{P{tm)P{tk)pP{tk)). (6.2) 

A sufficient condition for the validity of the semi-clasical approximation is therefore that the 
object 

Ak,m := TT{P{tjP{tk)pP{tk)) (6.3) 
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is much smaller than 1 for all tm > tk- 

Now assume for a moment that the semiclassical approximation holds, so that the class 



operators are given by Eq.(5.12). The probability of first crossing between tk and tfc+i is 
given by 

Tr(CfcpCi) = Tr(Cfcp) + 2Tr(pP(4+i))) - TT{pP{h)P{h+i))) - Tr(P(4+i)P(4)p). 

However since we are assuming that the semi-classical approximation holds, the last three 
terms cancel and we obtain 

TT{apcl) = TT{ap) = TT{pcl). 

Furthermore, consider a general off-diagonal term in the decoherence functional, where with- 
out loss of generality we take tk+i < tm 

Tr{CkpCl) = TT{P{tm)P{tk)p) - TT{P{tm)P{tk+i)p) 

-Tr(P(Wi)P(4)p) + Tr(P(t™+i)P(tfc+i)p). (6.4) 

We see that this vanishes if the semi-classical approximation holds. 

What we have shown therefore, is that the semi-classical condition and the decoherence 



conditions are closely related. If the object given in Eq.(6.3) is much smaller than 1 then 



both conditions are satisfied. This means firstly that the class operators are approximately 



given by Eq.(5.12), and secondly that the histories described by these class operators are 
approximately decoherent. 

In general, the probabilities for time of arrival as computed by decoherent histories differ 
from the heuristic ones obtained by the standard quantum mechanical analysis. We don't 



expect agreement since, as discussed in the Introduction, the heuristic formula, Eq.(l.l), is 
not in general of the canonical form required for genuine quantum mechanical probabilities. 
It is also the case that decoherent histories only ascribes probabilities to certain sets of 
histories. However if the object „ is small, then we have shown that the decoherent 
histories analysis reproduces the probabilities computed in the standard, heuristic, way. 

As well as clarifying the analysis of [9], note that all of the statements above apply to 
the case of a particle coupled to an environment. In this case the trace is to be taken over 
the system and environment, and P{tk) = P{tk)s ® le is a projection onto the degrees of 
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freedom of the system, s, tensored with the identity operator on the environmental degrees 
of freedom, e. This is the case of interest in this paper. 

Since the object A^^m plays a central role in our analysis it is interesting to ask if it has 
any physical interpretation. The answer is that it does. Recall that our class operator for 



crossing between tk-i and tk is given by Eq.(5.12). We see however that 



Pitk-i) - P(tfc) = P(tfe)P(4-i) - P(4)P(tfe_i). (6.5) 



The approximation leading from Eq.(5.11) to Eq.(5.12) is equivalent to dropping the second 



term on the right hand side of Eq.(6.5). The error in this approximation can be estimated 



by computing the probability associated with the class operator P(tfc)P(tfc-i). We see 
Tr([P(4)P(4_i)]p[P(tfc)P(4_i)]t) = Tr(P(tfc)P(4_i)pP(tfc_i)) = Ak-i,k. 



However Eq.(6.5) has a simple physical interpretation: it is the decomposition of the total 
current into right and left moving parts. The object A^^^ is therefore just the probability 
associated with the right moving current. Classically this is small by construction, since 
we will choose to work with wavepackets tightly peaked in negative momentum. Quantum 
mechanically however this term need not be small, indeed the existence of the backfiow 
effect [Ij shows that this term can sometimes be larger than the left moving current, even 
for wavepackets composed entirely of negative momenta. 

The semi-classical condition we are imposing, and by extension the decoherence condition, 
is therefore stronger than the condition for the current to be positive. Decoherence requires 
the absence of interference between crossings at different times, whilst the standard heuristic 
analysis includes some interference effects, provided they are not so large as to render the 
arrival time probabilities negative. The significance of this will be explored elsewhere. 

VII. DECOHERENCE OF ARRIVAL TIMES IN QUANTUM BROWNIAN MO- 
TION 

A. General case 

In the previous section we have shown that an arrival time distribution may be derived 
from decoherent histories, and that it agrees with the current provided there is decoherence. 
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We therefore turn now to the question of determining for which states and intervals decoher- 
ence is achieved. Recall that our aim is to show that the inclusion of an environment gives 
rise to decoherence of arrival time probabilities for generic initial states of the system. We 
will start our decoherent histories analysis with a discussion of the general case considered in 
Ref. [T5], but most of our detailed results will concern the near-deterministic limit discussed 
in the introduction, and in the discussion of the Wigner function above. 
The decoherence functional may be written in path integral form as 

V{a,a') = j^Vxj^ Vyexp (^^S[x] - ^S[y] + iW[x,y]^ Po{x,y), (7.1) 

where a, a' represent the restriction to paths that are, for example, in a; > at ti and in 
X < at ^2- Vr[a;,?/] is the Feynman-Veron influence functional phase which summarises the 
effect of the environment, and is given in the case of negligible dissipation by 

W[x,y] = '-^ I dt{x-y)\ (7.2) 

It is this phase which is responsible for the suppression of interference between paths x{t) 
and y{t) that differ greatly, and produces decoherence. 

The decoherent histories analysis we are about to perform was first attempted by Halliwell 
and Zafiris in Ref. [15]. Their conclusions are reasonable, but the analysis leading to them 
in fact contains a small error. This error arose due to a lack of appreciation of the role of 
the Zeno effect, as discussed in the Introduction. We aim to show here how the analysis 
may be modified in line with the treatment of the arrival time problem presented here and 
in Ref. [9]. 

We start from the expression for the density matrix for states that do not cross the origin 
in a time interval [0,t] (Ref.[l5], Eq.(4.34)) 



Prr{xf,yf) = Vx Vyexp 



^ f dt{x^ -f)-^ I dt{x - yf 



PQ{xo,yo). (7.3) 



Here the restriction is interpreted as x, ?/ > at times e, 2e.... The probability that the state 
does not cross the origin in this interval is then given by the trace of this expression. There 
are two interesting limits to be taken here. The first is letting e — )■ 0, and this recovers the 
notion of the restricted propagator. The second limit is that of strong decoherence effects, 
so that we can assume that the path integral is tightly peaked in {x — y), this recovers the 
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classical limit. The claim in Ref.[T5] is that we can take these limits simultaneously. To see 
why this is problematic we write the restricted path integral as a product of propagators, 



POO 

Prrixf,yf) = Y\ dxkdykJ{xk+i,yk+i,tk+i\xk,yk,tk)poixo,yo), (7.4) 
where tk = ke, to = 0, t„+i = t. Changing variables to X = x + y , ^ = ^{x — y) gives 
Prrixf,yf) = n / ^^i" / d^k(^)exp 



k=0 • 

De 
3h?' 



Po(Xo,eo). (7.5) 



Now in order for the product of propagators to be equal to a restricted propagator, we need 
to take e — > 0. However to recover the classical result we need to replace the limits of 
integration for ^k with ±oo, which requires that De » h'^/l'^ where / is some length scale. 
Clearly these two limits are incompatible for finite D. 

The standard discussion of the Zeno effect emphasizes the role of the so called "Zeno 
time", Tz = h{AH)^^. It is generally held that projections separated by times greater than 
this will not give rise to the Zeno effect. The issue is not without subtlety in the present 
case however, firstly because the Zeno effect is usually discussed in the context of unitary 
evolution, and this is not the case here, and secondly because characteristic times associated 
with the Zeno effect usually refer to selective measurement, and a projection onto positive x 
is certainly non-selective. This casts doubt on the role of the Zeno time as a timescale in this 
problem. Studies of the arrival time problem defined using a complex potential [7j suggest 
that the relevant timescale might be given by h/E. This issue will be taken up elsewhere 

[2H]. 

In any case it is clear that there exists a short timescale on which we can still approximate 
the limits in the ^ integrals above in the manner indicted, but that is short enough to give 
non-trivial crossing probabilities. Following the steps in Ref. ^L5\ then gives 

p{h,h) = l-TT{prr) = r dt f dp I dq^-^6{q)W[{p,q), (7.6) 

Jti J —oo J —oo 



m 



where the Wigner function, obeys, 

W[{p,q)= / dpodqoKr{p,q,t\po,qo,0)Wo{po,qo), (7.7) 
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where Kr is the "restricted" propagator, defined by 

p poo 

Kr{p,q,t\po,qo,0) = T\ dpk dqkK{pk,qk,tk\pk-i,qk^i,tk-i), (7.8) 
k=i 

with tk = ke. However because we are now dealing with an essentially classical system, this 
propagator is well approximated by its continuum limit. Understood in this sense, we see 
that the conclusions of Ref. [15] are in fact correct, even if the argument leading to them is 
not. The interesting subtleties that arise in this argument relate to the classical limit of the 
Zeno effect, and this will be explored elsewhere [25] . 

As important as this general case is, it is of interest to examine the simpler case of near- 
deterministic evolution, where the classical limit of the arrival time probability is simply 
expected to be the time integral of the current density. The analysis in this case simplifies 
considerably, and we will explicitly exhibit decoherence of histories for suitable intervals. 



B. The near-deterministic limit 



We have previously seen that, in the near-deterministic limit, if there is decoherence then 
the arrival time probabilities derived from decoherent histories agree with those computed 
from the current. We therefore turn now to discussing the conditions under which we have 
decoherence of crossing probabilities. The general picture we have in mind is illustrated 
in Fig.(l), we have an initial wavepacket defined at t = 0, evolved in the presence of an 
environment until time ti, and then we wish to compute the probability of crossing between 
ti and ^2- 



Our task is to compute the off-diagonal elements of Eq.(7.1), and this will be rather 
involved in general. However we are saved from having to do this by the observation in 
Section |Vl] about the relation between the semi-classical condition and the decoherence 



condition. It suffices therefore to compute the quantity Ak^m defined in Eq.(6.3), in the 
presence of an environment. We anticipate that this will be small simply from the form of 



Eq.(7.1). This is because the effect of the environment is to cause the density matrix to 
become tightly peaked around the classical path and classically, for p < the probability 
given by A^^^ is zero. We will see how this works in a specific example below. 
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*► X 



FIG. 1: The probability that a generic wavepacket arrives between ti and t2 can be expressed in 
terms of the state at ti or, via the propagator, in terms of an initial time t = 0. 

We take k = l,m = 2 without loss of generality, and we drop the subscript from now on 

A = Tr(P(t2)P(ti)pP(ti)) = jvxjvye^v (^^N - \s[y] -^j dt{x - yf^ Pt,ix^,y,), 

(7.9) 

where the histories a are those that start at Xi,yi < and finish at X2 > 0. In terms of the 



density matrix propagator, Eq.(|3.4|) 
A = 



I'D 

dx2 I dxi / dyiJ{x2,X2,t2\,xi,yi,ti)pt^{xi,yi) 



exp 



2h{t2-tij 

Transforming to new variables 



ITTl . . ,0 / \2\ ^(j^2 — ^1) / \2 \ / \ 

{{X2-X1} -{x2-yi)) [xi-yi) ] pt,{xi,yi). 



X= — ^ — , i = xi-yi, 



and writing the density matrix in terms of the Wigner function via 

p{x,y)= r dpe>(^~y^W{p,^y), 



(7.10) 



(7.11) 
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we obtain 



A 



m 



exp 



pX PCX} 

dX I dp 

X J-oc \2TTh{t2-ti^ 



(7.12) 



We see from this expression that there is a time scale, Ti = ^^2mh/ D, set by the environment 
on which locahsation effect are important. This timescale is the same as that on which the 
current becomes positive, as we saw earier. 

There are three cases to explore here, the first is where there is no environment, D = 0. 
This is the case covered in Ref.[3. The second case is the intermediate one, ti/r; >> 1 
but (^2 — << 1- This is the most general case in which we can expect to have 

environmentally induced decoherence. The final case is where ti/r/, {t2 — ti)/ri » 1. This 
is the case of very strong environmental coupling. 



For the free case, D = the integrals over ^ and X2 in Eq.(7.12) may be carried out to 
give. 



A 



/oo 
dpWt,iX,p)f 
'OO 



X 



mX 



where 



/(«) 



dy 



siny 



+ p 



y 



(7.13) 



(7.14) 



This expression for A is now identical to Eq.(5.43) of Ref.[7], and the same conclusions 
apply. We briefly repeat the analysis here however, since it turns out to be relevant for the 
second case. 

We note firstly that /(+oo) = 0, /(— oo) = 1, /(O) = 1/2. Now we assume that our state 
at ti is of the form 



W^.,(X,p) = -exp — -J^iv-V.) n 



(7.15) 



and that o"^ is large, so the state is tightly peaked in momentum. We can therefore integrate 
out p, setting p = po to obtain. 



dX 



1 



X 



mX 



h \{t2-h 



(7.16) 



In [7] it was noted that this type of integral is dominated by values of X for which 



< X < 



4iPor 



(7.17) 
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providied that Po(^2 — ti)/2m >> 1, and in this range the exponential term is approximately 
constant. We can approximate this integral then by intergrating from to — 7r/4|po|, taking 
X = in the exponential term and approximating, 

1 u 



f{u) 



This gives 



A 



TT h 



exp 



+ 0(m3). 
{Xo+poti/mf 



« 1, 



(7.18) 



(7.19) 



2 8(t|po| V 

there will therefore be decoherence for gaussian wavepackets tightly peaked in momentum, 
provided their momentum is such that Po(^2 — ti)/2m >> ie -Eo(^2 — ^i) >> h. In [7j it 
was argued that this conclusion also holds for orthogonal superpositions of gaussians. 



Turning to the intermediate case, since (^2 — ^i)/t/ << 1 we can set -D = in Eq.(7.12) 



and we obtain Eq.(7.13) again. Now however, Wt^{X,P) is the initial state evolved with 



an environment, and since ti/r; >> 1 this should be significant. To proceed we write the 



Wigner function at time ti in terms of the initial state at t = using Eq.(3.6) 



i^[Xo,Po] 



A 



dX / dpir(X,p,ti|Xo,po,0)/ 



X 



mX 



+ p 



(7.20) 



(7.21) 



Because the Wigner function propagator is a gaussian the analysis is similar to the first case. 
Recall that we are assuming t << and so p^/Dti » 1, this means we can integrate out 
p, setting p = Pq. This gives us 



F[Xo,Po] 



dXN 



exp 



{X-X,-p,U/mf)]f 



X 



mX 



h \it2-h) 



+ Po 



where, N,a,(3 are defined in Eq.(3.8). Now, 

pl/Dh » 1 

implies that 



1 



1 



pI n?i3 \h 



(7.22) 



(7.23) 



(7.24) 



so that for ti > n, 1/pl « l/f3h'^. This means, again like the case of D = 0, that the 
exponential term is roughly constant compared to /. Since we have Po(^2 — ti)/2m >> h, 
we follow the same procedure as the D = case, arriving finally at 



F[Xo,Po 



'IpoI 



exp 



/3 



{XQ+pQti/m) 



(7.25) 
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The value of A now depends on the relationship between the width of the initial state a and 
/3. However we can obtain an upper bound by ignoring the effects of the exponential term 
in F, this gives 

A<#5,/f=l,/?^f=)«l. (7.26) 



8bol V 4 ley pih \ti 

Finally we have the case of strong system-environment coupling. Since {t2 — ti) / ti » 1 



the integral over ^ in Eq.(7.12) will be peaked around ^ = 0, and we can therefore extend 
the limits of integration to (— oo, oo). This integral and the one over X2 may then be carried 
out and we obtain 







A = / dX I c/p-Erfc 



1. 



3m2 f^^p{t,-U) 



WtAX^p), (7.27) 



D{t2-hf\ m 
where Erfc is the complementary error function [22] • For large positive values of the argu- 



,2 



ErfcM ^ ^ ( 1 + O ( - ) 1 . (7.28) 



ment we have that, 

Erfciwl ^ 

Therefore, since the Wigner function is peaked around p < 0, A will be very small, provided 
(^2 ~ ^i) << Ts- This gives us an upper bound on the time interval, [ti,t2]; rather than a 



lower one. The lower time scale is provided by the condition ^2 ~ >> 'Ti = \/2mh/ D. 
This lower time scale is compatible with the condition that the current be positive. 

Note however that this lower limit is state independent. There will be states for which 
arrival time probabilities decohere on much shorter time scales than this, for example the 
simple cases which decohere in the absence of any environment will continue to do so in the 
presence of an environment, at least until a time ~ r^. 

The key point is that whether or not one can assign arrival time probabilities in decoherent 
histories depends on the form of the state at the time it crosses the origin. Environmentally 
induced decoherence produces mixtures of gaussian states from generic initial ones, and 
thus after a short time arrival time probabilities can be defined whatever the initial state. 
Crucially however it is not necessary for the system to be monitored whilst it is crossing 
the origin. The smallest time interval, 5t, over which we can define a decoherent arrival 
time probability is therefore set by the energy of the system and not the details of the 
environment and we must have E5t » h. This is in agreement with Ref. [7j, and also with 
the results of earlier works, concerning the accuracy with which a quantum system may be 
used as a clock [30] . 
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In conclusion then, for a general state, decoherence of histories requires that we evolve 
for a time much greater than t; = ■\j2mh/D before the first crossing time. This is because 
this is the time scale on which quantum correlations disappear and our initial state begins to 
resemble a mixture of gaussian states. After this time, we may define arrival time probabil- 
ities to an accuracy 5t, provided only that E5t » h. States which start as gaussians may 
be assigned arrival time probabilities without this initial period of evolution. This is in line 
with the general result that some coarse-graining is always required to achieve a decoherent 
set of histories in quantum theory. 

VIII. SUMMARY AND DISCUSSION 

In this paper wc have been concerned with deriving an arrival time distribution for open 
quantum systems, and comparing this with the classical result. We began by discussing the 
generalisation of the current, which is the classical arrival time distribution, to open quantum 
systems and in particular to the case of quantum Brownian motion. We found that in general 
the inclusion of an environment leads to extra terms appearing in the expressions for the 
current, compared with the those valid in the unitary case. However we have shown that in 
the limit of negligible dissipation these correction terms may be dropped. We then explored 
the resulting arrival time distribution and showed that it is non-negative after a time of 
order r^, and that after this time it can be written as the trace of the density matrix times 
a POVM. 

We then turned to the question of deriving this arrival time distribution from the decoher- 
ent histories approach to quantum theory. We extended the decoherent histories analysis of 
the arrival time problem to the case of a particle coupled to an environment. As expected, the 
inclusion of an environment produces decoherence of arrival time probabilities for a generic 
initial state. There are, however, some limitations to the permitted class of histories. For a 
generic state arrival times can only be specified after an initial time t » ti. Even after this 
time arrival times cannot be specified with arbitrary precision, coarse graining over intervals 
5t » h/E is required to ensure decoherence. We showed that the decoherence condition 
is very closely related to a semi-classical approximation for the evolution of the state, and 
that both conditions are satisfied if the time between projections is sufficiently large. This 
is a specific case of a more general connection between decoherence and classical behavior. 
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Our approach has proceeded at two levels. At the heuristic level the simple generalisation 



of Eq.(1.2) to the case of a particle coupled to an environment, Eq.(4.9), is a positive arrival 
time distribution after a time of order r^. This can also be written as the expectation value 



of a POVM, Eq.(4.14), and thus has the same form as other probabilities in quantum theory. 
On a more fundamental level, these expressions can be derived from the decoherent histories 
approach to quantum theory, where they are seen to be valid for times much later than r;. 

Although the arrival time probabilities computed from decoherent histories agree with 
the heuristic ones when we have decoherence, their derivation in this way represents a 
significant advance in our understanding. The great difficulty with regarding the current as 
the arrival time distribution is the arbitrary way in which one accepts these "probabilities" 
when they are positive, but declines to do so when they are not. Because decoherence is an 
essential part of the histories formalism, this arbitrariness is replaced with a consistent set 
of rules governing when probabilities may or may not be assigned. Whilst this may be of no 
consequence in the setting of a laboratory, it may prove hugely important in the analysis of 
closed systems, in particular the study of quantum cosmology [3T]. 

Another very interesting result we have presented is that the current becomes strictly 
positive after a finite time whilst assignment of probabilities in decoherent histories is only 
possible asymptotically. In some ways this difference between the heuristic analysis and 
decoherent histories is to be expected. Indeed, the current represents a linearly positive 
history [32], and it is essentially the condition of linear positivity that we have proven 
holds after a time of order r^. It is known that linear positivity is a weaker condition than 
decoherence [331 El] , and thus it is not surprising that demanding decoherence leads to a 
stricter limit on the assignment of probabilities than the heuristic analysis. It would be 
interesting to examine what is gained in this context by imposing decoherence rather than 
linear positivity. 

There are similarities here with the relationship between the current and Kijowski's ar- 
rival time distribution [35J, which may be shown to agree in the classical limit, but not 
more generally. Indeed because it is manifestly positive, one might expect the arrival time 
distribution computed from decoherent histories to be more closely related to Kijowski's 
distribution on shorter time scales. These issues will be discussed elsewhere. 
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